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Abstract 



One dimensional nano-beam photonic crystal cavities fabricated in silicon dioxide are considered in both simulation 
and experiment. Quality factors of over 10 4 are found via simulation, while quality factors of over 5 x 10 3 are found 
in experiment, for cavities with mode volumes of 2.0(A/n) 3 and in the visible wavelength range (600-716nm). The 
dependences of the cavity quality factor and mode volume for different design parameters are also considered. 



Silicon complementary metal-oxide-semiconductor 
(CMOS) compatible materials can be integrated with elec- 
tronics, and thus can have a large impact in computation, 
communication, and sensors. One important component of 
an optical system is a photonic cavity, able to store light and 
spectrally filter signals. The photonic crystal (PC) cavity is 
one type of optical cavity that allows for both high quality 
(Q) factors and low mode volumes {V mo d e ), having ubiq- 
uitous use in cavity quantum electrodynamics (cQED) [1], 
low threshold lasers [2], and optical control [3]. One of 
the most developed PC systems is the system of two di- 
mensional (2D) waveguides and cavities in a free-standing 
membrane, relying on total internal reflection for confine- 
ment in the third dimension. Silicon is transparent in the 
near infrared regime, and is heavily used for telecommuni- 
cation purposes near the wavelength of 1.5/xrri. However, 
silicon absorbs heavily in the visible wavelength range, and 
would be difficult to employ in light emitting and waveg- 
uiding devices at the visible wavelengths. 

Silicon dioxide (Si02, or silica), on the other hand, is 
transparent at the visible wavelengths, and similar to sili- 
con, is a promising material due to its low cost, compatibil- 
ity with electronics, and established fabrication techniques. 
The main hurdle in manufacturing photonic crystals in silica 
is its low index of refraction (n = 1.46). However, pertur- 
bation cavities that modulate the index of refraction in a 
waveguide system can be made even in low index materials. 
This approach has been used in creating high-Q cavities in 
2D PCs in both silicon [4-6] and low index materials such as 
silicon nitride (S13N4) [7], which has index n = 2.0, higher 
than that of Si02. Whereas a full photonic band gap in 
a 2D photonic crystal is difficult to achieve in low index 
materials, one dimensional (ID) nano-beam cavities that 
employ a gentle confinement technique enable high confine- 
ment with a small or incomplete band gap, while relying 
on total internal reflection in both directions perpendicu- 
lar to the beam length. Recent developments in ID nano- 
beam cavities with "potential well" designs have achieved 
the same Q-factors in silicon as in 2D photonic crystal cav- 
ities with comparable mode volumes [8,9], while also open- 
ing the door for high-Q cavities in SiaN4 for applications of 
nano-optomechanical coupled resonators [10] and coupling 



to active materials [11]. 

In this letter, we design and fabricate ID nano-beam 
cavities in silica. We follow the ladder cavity design used in 
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Figure 1: (a) The fabricated ID nano-beam cavity. Marker 
denotes lfim. (b) The electric field intensity (\E^ ) of the fun- 
damental mode supported by the cavity, (c) Band diagram for a 
beam with lattice constant a, w — 3a, d = 0.9a, h x — 0.5a, and 
h y = 0.7w, and another beam with the same parameters except 
for lattice constant a' = 0.9a. The dashed line indicates the light 
line in free space, (d) Design of the cavity, The plot shows the 
period (a) along the length of the beam as a function of N, the 
layer number counter from the center of the cavity. 

Si and S13N4 [10,12], but consider a silica slab suspended 
in free space, with lattice constant a, width w, slab thick- 
ness d, hole width h x , and hole height h y , as shown in Fig. 
QJb). We first obtain the band diagram of a periodic (or 
unperturbed) nano-beam waveguide using the three dimen- 
sional (3D) finite difference time domain (FDTD) method 
with Bloch boundary conditions. A sample band diagram 
is shown in Fig. |Tfc) for a beam with parameters: lattice 
constant a, w = 3a, d = 0.9a, h x = 0.5a, and h y = 0.7w, 
and also for a beam with the same parameters, except with 
lattice constant a' = 0.9a. As expected, the structure with 



the smaller lattice constant has slightly higher band fre- 
quencies, as this structure supports modes that have higher 
overlap with air. Because the lowest band edge mode of the 
structure with lattice constant a' (mode at the n/a' point) 
lies in the band gap of the structure with lattice constant 
a, it can serve as the defect mode in a beam with lattice 
constant a, which acts as the photonic crystal mirror. 
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where u> is the frequency of the mode, U is the total energy 
of the mode, and Pi is the power radiated in the i direction. 
We separate Q to t = l/(l/Q± + l/Q\\ + 1/Qpc), where Q± 
corresponds to radiation leaked in the z direction in Fig. 
QJb), Q|| corresponds to radiation leaked in the y direction, 
and Qpc corresponds to radation leaked in the x direction 
through the photonic crystal mirrors on the silica beam, 
and plot the results in Fig. ^a). We see that the gentle 
confinement method enables high reflectivity mirrors even 
in silica, as Q pc continuously increases with the addition of 
more PC mirror layers. In the case of a ID nano-beam, Q is 
limited by loss in the directions where the mode is confined 
by total internal reflection, namely Q± and Qn. As seen 
in Fig. EJJa), the limiting factor in Qtot is Q± in this case. 
Thus, Qtot could be increased by improving the design of 
the periods that correspond to the photonic crystal cavity, 
possibly by parameter search, genetic algorithms [14], or 
inverse designs [15]. 



Figure 2: (a) The direction specific Q-f actors of the nano-beam 
cavity as a function of the number of photonic crystal mirror lay- 
ers surrounding the cavity. With respect to Fig. \j](b), Q± corre- 
sponds to radiation leaked in the z direction, corresponds to 
radiation leaked in the y direction, and Qpc corresponds to ra- 
diation leaked in the x direction. Qtot is the parallel sum of Q±, 
Qn, and Qpc- (b) The Q -factors and (c) the mode volumes of 
cavities with the same air hole design, but different beam widths 
and thicknesses. The reference dot sizes are for Q — 2.0 x 10 4 
and V mo de = 2.0(A/n) 3 in (b) and (c), respectively. 

Next, we use the perturbation design suggested by pre- 
vious references and introduce a parabolic relationship be- 
tween the lattice constant and the z-coordinate, thus form- 
ing an optical potential well [13]. In particular, we choose 
a minimum effective lattice constant of 0.92a at the center 
of the potential well, and the remainder of the lattice con- 
stants as shown in Fig. [T^cl) . The perturbations of lattice 
constants span 7 periods away from the center of the cavity. 
We simulate the full cavity structure again using the FDTD 
method, and obtain the fundamental TE-like mode with 
frequency a/A = 0.454, Q = 1.6 x 10 4 , V mode = 2.0(A/n) 3 , 
and electric field intensity (|-E| 2 ) shown in Fig. QJb). This 
represents a more than one order of magnitude increase in 
Q-factor and a 7 fold reduction in mode volume compared 
to a 5/im diameter silica microdisk cavity with the same 
silica thickness. 

We also simulate the effect of the number of photonic 
crystal mirror layers on the cavity Q. In particular, we 
define the quality factor in the direction i as Qi = cull/ Pi, 
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Figure 3: (a) The cavities spectra measured in reflectivity from 
structures with different lattice constants are normalized and 
shown together, with fits to Lorentizian lineshapes plotted on top 
of the data points. The cavity spectra are shown from left to 
right with increasing a. (b) The Qs of the cavities shown in part 
(a) plotted against the wavelengths of the cavities, (c) The an- 
gle dependence of the reflectivity amplitude. The horizontal axis 
corresponds to the half-waveplate angle, which is placed in front 
of the objective lens that is in front of the chip ( and thus placed 
in both the incident and collection paths). The fit to the reflec- 
tivity amplitude shows a period very close to tt/4, indicating a 
linearly polarized cavity mode. 

Furthermore, we simulate the cavity with the same pat- 
tern of holes in the x-dircction, while changing the beam 
width w and thickness d, keeping h x = 0.5a and h y = 0.7w. 
The resulting limiting Q (the parallel sum of Q± and Qu) 
and Vmode are shown in Fig. HJb) and (c), respectively. 
We see that the both the Q-factor and V mo de increase with 
the slab width and thickness, which is expected as larger 
cavities have higher confinement but higher mode volumes. 
In order to employ this cavity in cQED applications, we 
need to maximize the Q/V mo de ratio. Such a maximum is 
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achieved with w = 2.6a and d = 1.1a, with Q = 2.0 x 10 4 
and V mode = 1.8(A/n) 3 . 

We next fabricate the cavity with the same dimensions 
as above. We start by oxidizing a silicon wafer, forming a 
270nm layer of SiC>2 on top of silicon. Next, e-beam lithog- 
raphy is performed with a 250nm layer of ZEP-520A as the 
resist. After development of the resist, the pattern is trans- 
ferred to the oxide layer with a CHTVC^ chemistry dry 
etch. Finally, the beam is undercut with a XeF2 etcher, re- 
moving approximately 4/im of silicon under the oxide layer. 
The final fabricated structure is shown in Fig. [Ha). We 
vary the lattice constant a in fabrication to create cavities 
with a variety of wavelengths. 

We characterize the cavities using the cross-polarized 
reflectivity measurement technique [16, 17]. In summary, 
white light linearly polarized 45° from the cavity polariza- 
tion is directed at the cavity through an objective lens from 
above (z direction), and emission in the polarization orthog- 
onal to the excitation (135°) is collected also from above 
and detected by a spectrometer. The fundamental mode of 
the cavity is linearly polarized in the direction perpendicu- 
lar to the beam length (y-polarized in FigQJb)), much like 
the fundamental TE mode of a rectangular slab waveguide. 
We observe cavities of different resonant wavelengths that 
span 600-716nm in the visible wavelength range for differ- 
ent lattice constants, as shown in Fig. [3](a). The cavities 
are shown from left to right with increasing a. The mea- 
sured cavity wavelengths and Qs of the cavities are plotted 
in Fig. [3£b). We also observe a slight decrease in Q with 
increasing a, which is expected from the FDTD simulated 
trend of decreasing Q with decreasing d/a. By placing a 
half-waveplate in the incident/collection path, we confirm 
that the observed cavity mode is linearly polarized, as the 
angle dependent reflectivity amplitude has a period of 7r/4 
with respect to the half-waveplate angle (Fig. OJc)). 

In conclusion, we have demonstrated a photonic crystal 
cavity in silica with theoretical Q-factors above 10 4 , and ex- 
perimental Q-factors above 5 x 10 3 , while maintaining low 
mode volumes of 2.0(A/n) 3 . Such a system could be en- 
hanced by improving the cavity design for higher Q-factors 
through manipulation of the cavity hole placement. In addi- 
tion, the coupling of emitters to this low index system could 
be explored. Finally, the silica cavity could be employed to 
couple optical fields and mechanical oscillations. 
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